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ABSTRACT 

A mathematical model of a hydrogen/oxygen alkaline fuel cell is presented that 
can be used to predict the polarization behavior under various power loads. The major 
limitations to achieving high power densities are indicated and methods to increase 
the maximum attainable power density are suggested. These performance indications 
can help future research and the design of alkaline fuel cells. 

The alkaline fuel cell model describes the phenomena occurring in the solid, liq- 
uid, and gaseous phases of the anode, separator, and cathode regions based on porous 
electrode theory applied to three phases. Fundamental equations of chemical engi- 
neering that describe conservation of mass and charge, species transport, and kinetic 
phenomena are used to develop the model by treating all phases as a homogeneous 
continuum. Gas phase diffusional resistances are considered by calculating the spatial 
variation of the partial pressures of oxygen, hydrogen, and water vapor in the gas 
phase. The liquid phase diffusional resistances are accounted for by considering the 
concentration distributions of dissolved oxygen and hydrogen in KOH. The variation 
of the KOH electrolyte concentration is also accounted for by including the ionic 
resistance effects. Electronic resistances are considered by calculating the solid elec- 
trode potential drops in the porous gas diffusion electrodes. By developing a complete 
model of the alkaline fuel cell, the interaction of these various resistances can be 
investigated under conditions that simulate actual fuel cells. 

A sensitivity analysis of the various transport and electrokinetic parameters indi- 
cates which parameters have the most influence on the predicted current density and 
over which range of potentials these parameters affect the fuel cell performance the 
most. This information can be used to decide which parameters should be optimized or 
determined more accurately through further modeling or experimental studies. The ef- 
fect of various design parameters on the limiting current density are also investigated 
to determine if optimal values exist for the parameters. These parameter sensitivi- 
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ties and optimal design parameters can help in the development of better three-phase 
electrodes and separators for the alkaline fuel cell. 
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CHAPTER I 
INTRODUCTION 

The alkaline fuel cell (AFC) is capable of providing a clean, efficient, and high 
powered source of electrical energy. The relative ease of operation, low weight and 
volume, and reliable performance has made the alkaline fuel cell an attractive power 
source especially for the American space program as well as for electric vehicles, 
defense, stationary power sources, portable generators, and small submersibles to name 
a few (1-3). However, for applications that require even higher power densities, the 
performance of the alkaline fuel cell needs to be improved. A mathematical model 
of the alkaline fuel cell can assist in understanding better the phenomena occurring in 
the system as well as help in the design of fuel cells. 

A. Fuel Cells 

A fuel cell is a reactor where a fuel and an oxidizer react electrochemically 
producing products and releasing energy. Typically, during a combustion process, the 
energy of the fuel is released as heat which is then converted to electrical energy 
through a generator. However, one of the main advantages of the fuel cell is that 
the chemical energy of the fuel is directly converted into electrical energy which 
can be used to produce work. By avoiding the chemical to mechanical conversions, 
efficiencies of 50 to 60% can be achieved in fuel cells in contrast to the Carnot limited 
efficiencies of 15 to 25% in combustion engines (2). 

Recent research in fuel cells has led to the development of various types of fuel 
cells such as the phosphoric acid, molten carbonate, solid oxide, and solid polymer 
electrolyte fuel cells as discussed by (1,4—7). However, these systems still cannot 
compete with the alkaline fuel cell’s high power density. Continued research in 
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electrocatalysts, porous gas diffusion electrodes, and materials has increased the power 
density for the AFC dramatically since the 1960’s (8). By developing a mathematical 
model of the alkaline fuel cell, the major factors that limit the performance of the 
system can be investigated. These performance indications can be used to guide 
future research in alkaline fuel cells. 

B. Objectives 

The maximum power density obtained from alkaline fuel cells has undergone many 
advances in recent years due to improved catalysts and electrode materials and also 
due to optimized operating conditions and fuel cell design. It has been known that in 
order to increase the power density in AFCs, the activation, concentration, and ohmic 
polarization should be minimized (4). For example, improved electrocatalysts for the 
oxygen reduction reaction in alkaline electrolytes has helped reduce the activation 
polarization (9); increased solid and solution conductivity has helped reduce the 
ohmic polarization (10); dual porosity electrodes have helped reduce concentration 
polarization (10-12). However, in order to investigate these phenomena, numerous and 
expensive experimental tests need to be performed. Furthermore, the interaction among 
the three types of polarization may cause difficulties in isolating the characteristics of 
a particular type of polarization. Mathematical modeling can help determine how 
changes in parameters and operating conditions will influence the various types of 
polarization which subsequently affect the performance of the fuel cell. The AFC 
model can help identify parameters and concepts that limit the performance of the 
fuel cell based on today’s state-of-the-art technology. Additionally, the model can be 
used to investigate the effects of hypothetical advances in technology on the predicted 
performance. 

The objectives of this work are to develop a realistic mathematical model of a 
hydrogen/oxygen alkaline fuel cell and to use the model to predict the maximum 
attainable power density. Fundamental equations of chemical and electrochemical 
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engineering are used to describe conservation of mass and charge, species transport, 
and kinetic phenomena in the fuel cell. Gas and liquid phase diffusional resistances as 
well as ionic and electronic resistances in the solution and solid phases, respectively, 
are accounted for in the model. The production and removal of water are also 
considered so that all electrochemically produced water will leave the system with 
the excess gas streams. The effects of the various transport, thermodynamic, and 
kinetic parameters on the fuel cell’s performance are investigated to identify the 
more influential parameters. The mathematical model is then used to optimize the 
thicknesses and porosities of the fuel cell’s various regions by maximizing the power 
density. 

C. Structure of Report 

This report contains the development and analysis of an alkaline fuel cell model. 
The motivation for developing this model is presented in Chapter I. A review of the 
relevant literature on porous gas diffusion electrodes and alkaline fuel cells is discussed 
in Chapter II. The development of the alkaline fuel cell model is shown in Chapter 
III along with some results of the dependent variables and model predictions. The 
sensitivity of the model predictions and the optimization of certain features of the fuel 
cell are presented in Chapter IV. The results of the alkaline fuel cell model are then 
summarized and recommendations for future studies are given in Chapter V. 
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CHAPTER II 

LITERATURE SURVEY 


A. Three Phase Electrodes 

The need to increase current densities in electrochemical systems has led to the de- 
velopment of three phase electrodes (porous gas-diffusion electrodes) containing large 
interfacial surface areas between the solid electrocatalysts, electrolyte, and gaseous 
pores. However, since these electrodes contain a tortuous and non-uniform distribu- 
tion of catalysts, gaseous-filled pores, and liquid-filled pores, it is difficult to accu- 
rately describe the phenomena occurring in the electrodes. Numerous experimental 
and modeling studies have been conducted over the years to describe and enhance 
the performance of these porous electrodes as discussed by (13-17). The major dif- 
ficulty in modeling three-phase porous electrodes is in describing the electrode struc- 
ture. Most gas-diffusion electrodes incorporate a hydrophobic agent such as Teflon or 
polytetrafluoroethylene (PTFE) mixed with the electrochemically active hydrophilic 
catalysts. This creates a stable three-phase boundary but also creates a very complex 
wetting phenomena making it difficult to characterize the structure of the electrode. 

One of the earliest models for a three phase electrode was developed by Will 
(18,19) who conceptualized the three phase electrode as consisting of a thin liquid film 
in the form of a meniscus covering the electrocatalyst. Will concluded that most of the 
current is generated in a small diffusion region near the thin film/meniscus boundary. 
A slight variation to this approach was shown by Rockett and Brown (20) who allowed 
the thin film to have a variable thickness. These models only included solution phase 
diffusional and ohmic resistances, neglecting any gaseous and electronic resistances. 
The interaction of the KOH and dissolved O 2 kinetic parameters were shown to be 
responsible for the limiting current density in Rockett’s model, rather than just a single 
parameter as concluded in prior models. A more complete thin film/meniscus model 
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was developed by Bennion and Tobias (21) where they considered the diffusion and 
migration of the relevant species and the solubility of O 2 in the electrolyte. Their 
model results showed that the current density is controlled by the charge transfer 
overpotential and the ohmic resistance drop in the thin film. They also concluded that 
the transport of dissolved O 2 is not rate limiting in the thin film meniscus region but 
is limiting in the bulk electrolyte regions. 

Iliev et al. (22) developed a simple model of a gas-diffusion electrode that 
considered the Knudsen diffusivity of gaseous O 2 , gas dissolution, and ohmic drop 
only. They considered the diffusion of gaseous O 2 as occurring by Knudsen diffusion 
rather than molecular diffusion since the pore radii were on the order of 30 nm for 
their work. They concluded that the porous structure of the electrode and the mode 
of mass transfer through the gaseous pores are responsible for the electrochemical 
activity in the electrode. Unfortunately, their model did not include ionic or electronic 
resistances which can drastically influence the system. 

The models by Will, Rockett, and Bennion are usually termed microscopic mod- 
els since they try to describe the physical phenomena occurring at the three-phase 
boundaries. Another approach to describing three-phase electrodes is to treat them 
macroscopically where the unknown structural geometries of the porous electrodes are 
assumed to be homogeneously mixed. One such approach is the agglomerate model 
(23—25) where the catalyst particles and the gaseous and liquid filled pores co-exist in 
a homogeneous continuum. One of the earlier agglomerate models was developed by 
Giner(23) who accounted for the diffusion of a dissolved gaseous species in parallel 
with an electrochemical reaction in a cylindrical agglomerate. Giner also accounted 
for the solution phase potential variation. The model results showed the performance 
of the electrode as a function of agglomerate radius and catalyst utilization. 

Cutlip (26) developed an analytical model of the mass transfer process occurring in 
gas diffusion electrodes. This model was primarily designed to study the mass transfer 
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effects due to low concentrations taking into account gaseous diffusion, gas dissolution, 
and the transport of dissolved gas. Catalyst effectiveness factors (27) were used to 
modify the electrocatalyst based on different shape factors. Cutlip concluded that the 
limiting current density is not influenced by the liquid phase diffusional resistances. 
Another important conclusion from Cutlip’s model was that the coupling of the gas- 
phase diffusional resistance and electrode thickness were strongly influential on the 
limiting current density. The limiting current density was found to be proportional 
to electrode thickness when the gas-phase diffusional resistance was non-existent. 
Conversely, the limiting current density was found to be independent of electrode 
thickness when gas phase diffusional resistances were present. 

An improvement to Giner’s flooded agglomerate model and to Cutlip’s earlier 
model was presented by Iczkowski and Cutlip (28) who developed a fairly complete 
model of a three-phase electrode. This model accounted for the gaseous diffusion 
resistances, solution phase diffusional resistances, and solid and liquid potential 
variations. They fitted their model to experimental polarization curves by adjusting the 
radius of the agglomerates, electrolyte film thickness, and the porosity-tortuosity factor 
of the gaseous and liquid filled pores. They applied their model to a phosphoric acid 
electrode up to current densities of 370 mA/cm 2 and concluded that ohmic resistance 
was the major cause of polarization losses in the electrode followed by Knudsen and 
molecular diffusion effects. However, they concluded that the solution phase resistance 
of dissolved oxygen diffusion contributed the least to the polarization losses. 

An agglomerate model for a double-layered oxygen electrode was also given by 
(29) where it was assumed that the electrochemical reaction was analogous to the 
diffusion and reaction of a gaseous reactant in a porous catalyst pellet surrounded 
by a gaseous film. Unfortunately, too many simplifying assumptions ( e.g ., constant 
overvoltage, constant electrolyte concentration, simple Tafel kinetics, etc.) limit the 
applicability of this approach to more complex phenomena occurring in the electrode. 
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A different macroscopic model applied to two-phase porous electrodes was devel- 
oped by Newman and Tobias (30) which is commonly referred to as porous electrode 
theory. This approach treats the different phases in the electrode as a homogeneous 
continuum that can be characterized by measurable quantities such as the porosity and 
specific surface area as described by (17). This macroscopic approach to describing 
the porous electrode assumes that the variables of interest are continuous in time and 
space allowing them to be averaged over a small volume element in the electrode. A 
macroscopic treatment of the chemical, electrochemical, and physical processes occur- 
ring in a three-phase electrode was developed by Darby (31) and extended by White 
et al. (32) based on porous electrode theory. It is, perhaps, rather unfortunate that 
the continuum approach has been criticized as being too abstract in comparison to 
the agglomerate approach (23). A close examination of the two macroscopic methods 
show that both methods are essentially describing the same type of conceptualization; 
the major differences being the way electrode surface areas and film thicknesses are 
handled. 

B. Alkaline Fuel Cell Models 

As shown in the previous section, there have been many investigations into three 
phase electrodes, in particular, the oxygen electrode since it is responsible for most 
of the polarization losses in AFCs. Unfortunately, no complete model of the alkaline 
fuel cell exists in the open literature that contains both electrodes and the separator. It 
is difficult to generate any firm conclusions on an entire fuel cell assembly when only 
a single electrode model is used since the interactions between the anode, cathode, 
and separator are not considered. 

All of the single electrode models presented above simplified or neglected various 
forms of resistance ( i.e gas phase diffusion, liquid phase diffusion, electronic, or ionic 
resistances). Some of these simplifications such as constant electrolyte concentration 
or constant solid potential are, as will be shown later, justified when a low overvoltage 
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is applied as was done in these models. However, at higher overvoltages, the 
assumptions of constant potentials and concentrations fail since large polarization 
effects are occurring. Hence, the previous models are not sufficient when trying 
to predict high power density performance. By considering all pertinent regions in the 
alkaline fuel cell, performance predictions can be made at high power densities. 
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CHAPTER III 

A ONE-DIMENSIONAL MODEL OF AN ALKALINE FUEL CELL 

A. Introduction 

A schematic of an overall alkaline fuel cell system with its associated gas channels 
is shown in Fig. 1 which is used as a basis for the model development. The modeling 
region itself is presented in Fig. 2 showing a conceptualization of the three-phase 
electrodes. The AFC operates by flowing either dry or humidified hydrogen and 
oxygen gases through the anode and cathode gas flow channels, respectively. As 
these gases flow through the channels, humidified hydrogen gas diffuses into the 
gas diffusion region of the anode while humidified oxygen gas diffuses into the 
gas diffusion layer of the cathode. Next, the gases diffuse into the reaction layers 
of the electrodes, where the gases further diffuse in the gaseous phase as well as 
dissolve into the KOH electrolyte. In the middle of the fuel cell is a non-conducting 
microporous separator matrix which is assumed to prevent any gases from diffusing 
across the system. In some alkaline fuel cells, the electrolyte is circulated out of the 
system which helps maintain a constant electrolyte concentration and assists in heat 
and water removal. A design constraint for this work is that the KOH electrolyte 
does not circulate outside the system. This will cause large electrolyte concentration 
variations across the electrodes and separator. In practical fuel cells, the liquid water 
produced in the system dilutes the KOH electrolyte to a certain extent. It is assumed 
that all electrochemically produced liquid water evaporates into the gas streams. This 
assumption is reasonable since in actual alkaline fuel cells the average electrolyte 
concentration reaches a constant value at steady state. 

The design of the porous gas-diffusion electrodes is critical in order to achieve high 
performance in the alkaline fuel cell. It is important that the three-phase boundaries in 
these electrodes provide a high utilization of the catalyst clusters which will give large 



10 



Cathode Gaa 
Channel 


Current 

Collector 

Cathode Gas 
Diffusion Layer 

Cathode Reaction 
Layer 


Separator 


Figure 1. Overall schematic of a hydrogen/oxygen alkaline fuel cell. 
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Figure 2. Detailed schematic of the three-phase 
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surface areas for the electrochemical reactions. It is also important that the electrodes 
have good electrical conductivity and are corrosion resistant. The preparation and 
design of some high performance gas-diffusion electrode are given by (10,33-38) 
where it has generally been concluded that the reaction layer should be optimized 
to contain a large number of catalytic clusters in the electrolyte and a large number 
of gas dissolving sites. Typically, multilayer electrodes are used as shown in Fig. 2 
where each electrode contains a gas diffusion region consisting of hydrophobic gaseous 
pores only, and a reaction region consisting of hydrophobic gas-filled pores as well as 
liquid-filled pores. The hydrophobicity in the electrodes is obtained by impregnating 
the porous electrodes with a wet-proofing agent such as teflon, polytetrafluoroethylene 
(PTFE), or wax which also serve as a binding material for the electrodes (10). The 
quantity and distribution of these wet-proofing agents have a strong effect on the fuel 
cell performance as shown by (11,39). 

As the dissolved gases diffuse through the electrolyte, they electrochemically react 
according to the following reactions: 

Cathode : 0> + 2 1 1 2 O ( / j + 4e~ — > 40H [1] 

and 

Anode : H -2 + 2011“ — » 2H20(q + 2e“ [2] 

Hence, the overall reaction is the production of water 

Total : H? + — O 2 — > [3] 

with the simultaneous liberation of electrical energy. There have been numerous in- 
vestigations into the mechanism of the oxygen reduction reaction in alkaline solutions. 
The two more common mechanisms for oxygen reduction involve the production of 
peroxide as either an intermediate species (40-43) or as a reaction product (44,45). 
However, the kinetic parameters vary depending on the type of electrode substrate and 
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electrocatalyst (9,46-49) as well as on operating conditions such as concentration, pH, 
temperature, and oxygen partial pressure (50,51). Therefore, since too many factors 
influence the oxygen reduction mechanism, it will be assumed that the direct four 
electron transfer process, as given by Eq. [1], occurs in the oxygen electrode. 

B. Phenomenological Equations 

In order to investigate the performance of an alkaline fuel cell, the phenomena 
occurring in the separator and in the three-phases of the electrodes need to be described. 
Including all regions of the cell allows the anode and cathode to interact through the 
continuous distribution of the solution phase potential and the electrolyte concentration 
across the system. This interaction is impossible to investigate with a single electrode 
model. To account for the complexities of the microporous structures, porous electrode 
theory (17), is used as a basis to describe the electrodes. This theory allows the 
superimposition of two or more phases into a single, homogeneous continuum. Thus, 
the gas diffusion layers are described by a homogeneous continuum of gaseous-filled 
pores and solid electrode material. Similarly, the reaction layers can be described 
by superimposing the gaseous-filled and liquid-filled pores with the solid electrode 
particles. 

To develop a complete model of the alkaline fuel cell, the various forms of 
resistance need to be considered. In the gas phase, the reactant gases (H 2 and O 2 ) 
diffuse via molecular diffusion through water vapor in their respective electrodes 
contributing to the gas phase diffusional resistances. Thus, the partial pressures of 
hydrogen (Ph 2 ), oxygen (Po 2 )> and water vapor in the anode (T^o) ar| d cathode 
(Pfj q), need to be determined. Electronic resistances occur through a potential drop 
in the anode ( E a ) and cathode (E c ) while an ohmic resistance arises due to a varying 
solution phase potential (</>). To account for the liquid phase diffusional resistances, 
the concentrations of dissolved hydrogen (Ch 2 ) and dissolved oxygen (Co 2 ) need 
to be determined. The electrolyte concentration ( C e ) varies in the solution phase 
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contributing to the ionic resistance. Additionally, the volume average velocity ( v ) 
needs to be determined in order to consider convective effects in the solution phase. 
The fundamental equations needed to solve these eleven dependent variables will now 
be described in their general form. After these generalized equations have been applied 
in each region of the fuel cell, the necessary boundary conditions will be developed. 

A one dimensional mathematical model of the alkaline fuel cell can be developed 
by considering conservation of mass and charge, transport of species, and reaction 
kinetics in each of the regions of the fuel cell as shown in Fig. 2. The equation of 
continuity for species i can be written in the general form for a porous medium 

r)fC 

^1= -V-N , + + (i = 0 2 ,H 2 ,H 2 0,+,- o) [4] 

where R v t and R e t represent rates of production for material that is “brought in” from 
across a phase boundary and for material that is produced from an electrochemical 
reaction, respectively. Since porous electrode theory is being used, the reaction rate 
terms are included throughout the electrode, as indicated in Eq. [4], instead of 
being treated as a boundary condition as done in agglomerate models. Note that 
the ,and o (water) species in Eq. [4] correspond to the K + ions, OH - ions, 
and solvent, respectively. Furthermore, the V operator in Eq. [4] and in subsequent 
equations is represented by ^ where z is in the horizontal direction in reference to 
Fig. 2. The flux expression for species i, N,, depends on whether species / is in the 
gas or solution phase. In the gas phase, the Stefan-Maxwell equation can be simplified 
for a binary gas mixture of species i and j (52): 

9 

- /?T 

Vy, = E - W N <) (M = 0 2 or H 2 , H 2 0) [5] 

j = l 1 U{ ’J 

where the effects of diffusion and convection are accounted for. The ionic flux 
expression (53) 

N, = -P'YC, - z lUl FC t V<j> + Civ {i = 0 2 , H 2 , o) [6] 
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can be used to represent the transport of species i in the liquid-filled pores of the porous 
electrodes. This expression accounts for the diffusive, migration, and convective 
effects in the solution phase through the first, second, and third terms, respectively. 
Note that Eqs. [5] and [6] contain effective diffusivities, V t , which are related to the 
free stream diffusivities by a porosity and tortuosity factor: 


r 


and the Nemst-Einstein relation is assumed to relate the mobility, to the diffusivity 
of species i. 


Ui - 


A 

RT 


[ 8 ] 


The electrochemical reaction rate per unit of electrode volume, is expressed 
for species i in the form: 


R\ = 


sni l i 

nF 


[9] 


where the stoichiometric coefficients, s t , are given by expressing the electrochemical 
reactions in the form 


V siM z ' — > ne [10] 


The local current density, /, is expressed by the Butler- Volmer electrochemical rate 
expression 
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Ci 

c: 


TTr exp 


oc a nF 

RT 
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Ci 
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exp 


-a c nF 

RT 


[ 11 ] 


where the overpotential, 7/, is given by 


7/ = E - <f> - U ref 


[ 12 ] 
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Note that Eq. [11] is used for the hydrogen oxidation and oxygen reduction reactions 
rather than Tafel expressions in order to account for the effects of the reactant and 
product concentrations over the entire range of potentials investigated. 

The rate of production of hydrogen and oxygen gas into the electrolyte across a 
phase boundary, R v { , is approximated by 


ll v = -a 9 V\ 

X 1 



[ 13 ] 


where //, is Henry’s law constant in mol/(cm 3 atm) for species i, 6 is the diffusion 
layer thickness, and a 8 is the specific surface area of the gaseous pores. Note that the 
term, a 9 V l J8, could be replaced by a mass transfer coefficient, but this may lessen 
the appeal of the functionality of the dissolution rate expression. This rate expression 
assumes that equilibrium will be established at the gas-liquid interface following a 
Henry’s law expression. 

These equations describe the physical phenomena believed to be occurring in the 
alkaline fuel cell and will now be applied to the specific regions of the fuel cell. 


Gas Diffusion Regions 

Each electrode contains a gas diffusion layer to prevent the electrolyte from 
weeping into the gas stream and to provide structural support to the electrode. It 
may be assumed that the hydrophobicity of the electrode will prevent any liquid 
from entering this region. Therefore, only hydrogen gas and water vapor will exist 
in the anode gas diffusion layer and only oxygen gas and water vapor will exist in 
the cathode gas diffusion layer. Additionally, the solid electrode material in the gas 
diffusion regions will experience an ohmic drop that follows Ohm’s law. Hence, 
each gas diffusion region has three unknown variables: P \\ 2 , anc * ^ or 

anode layer and Pq 2 , P£ 2 q, ar >d E c for the cathode gas diffusion layer. In the gas 
diffusion regions, the reactant gases (H 2 or O 2 ) diffuse through water vapor from the 
gas channel/gas diffusion interface to the gas diffusion/reaction interface. The water 
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vapor in these regions comes from two sources: water vapor that enters the system 
with the reactant gases and the water vapor that evaporates from the electrolyte in 
the gas reaction regions. Thus, the water vapor itself can diffuse in or out of the 
gas diffusion regions depending on the inlet and reaction conditions. To properly 
account for the relative fluxes of reactant gas to water vapor needed in the Stefan- 
Maxwell expression, Eq. [5], a water balance is needed about the system. From the 
stoichiometry of the overall reaction, Eq. [3], the flux of hydrogen gas is related to 
the total flux of water leaving the system 

N g H2 = -Nj h0 [14] 


and similarly for oxygen gas: 

N Ch = ~2 N h*° 

Since the total flux of water vapor is simply the sum of the water fluxes that leave 
with the anode and cathode gas streams, 


N 


T 

H 2 0 


= Nh,o + N 


h 2 o 


[16] 


then a water fraction, f a , can be defined that relates the amount of water that leaves 
through the anode to the total amount of water generated (54): 


/« = 


N'h 2 o 


^h 2 o + n h 2 o 


[17] 


Similarly, a cathode fraction (f c = 1 - f a ) can be defined which gives: 


fa = 1 - 


* e H t O 


* ko + N 


HiO 


[18] 


Equations [14], [16], and [17] can be combined to relate the flux of hydrogen gas to 
the flux of water vapor in the anode: 


N 


h 2 o “ 


-/.A'ff, 


[19] 



Similarly, Eqs. [15], [16], and [18] can be combined to give the relationship between 
the oxygen and water vapor flux in the cathode: 
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NjhO = -2(1 - f.Wo, 


[ 20 ] 


These flux ratios can be inserted into the convective term of the Stefan-Maxwell flux 
expression, Eq. [5], and combined with the equation of continuity, Eq. [4] (where R ? 
and R\ are zero since there are no reactions in this region), to give: 


dt 


= V • 


v e 2 p t 

P}] 2 0 + fa^lh 


VP,,. 


[ 21 ] 


in the anode and 


dt*Po 2 

dt 


K,Pt 


pc 

1 H 2 


o 


+ 2(1 - f a ) Pq 


-VPo 2 


[ 22 ] 


in the cathode. Since no production or consumption of any kind occurs in the gas 
diffusion layers, the total pressure will be constant in each region 


1 T 


P\l 2 + ^H 2 0 


[23] 


in the anode and 


P C T = Po 2 + Pr. 


o 


[24] 


in the cathode. 

The last dependent variables that need to be determined in the gas diffusion regions 
are the anode and cathode solid potentials, E a and E c , respectively. The ohmic drop 
in these regions can be described by Ohm’s law: 

VE=-- [25] 

a 

where a is the conductivity of the electrode and I is the total cell current density. 
Since there are no electrochemical reactions in the gas diffusion region, then the 
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current density is constant so that the gradient of the current density is zero. Hence, 
the ohmic drop in the anode gas diffusion layer is given by 


V 2 E a = 0 [26] 

and similarly in the cathode gas diffusion region: 

V 2 E c = 0 [27] 

The steady-state forms of the governing equations for the anode gas diffusion region 
are summarized in Table I where the spatial dimension was made dimensionless by 
setting: 

f - z 1281 

Similarly, the steady-state forms of the cathode gas diffusion region’s governing 
equations are shown in Table II. 

Gas Reaction Regions 

In addition to containing a gas diffusion layer, each electrode has a gas reaction 
layer where the electrochemical and dissolution reactions occur. The reactant gases 
(O 2 or H 2 ) diffuse through water vapor in the gas pores of the electrode while some 
of the gases dissolve into the liquid-filled pores. The dissolved gases further diffuse 
in the solution until they reach a reaction site where they react electrochemically. 
The electrochemical reactions are influenced by the electrolyte concentration, soludon 
potential, solid electrode potentials, volume average velocity, and dissolved reactant 
gas concentrations. Thus, in the anode gas reaction region there are seven variables 
to solve for: C e , <f>, E a , C\\ 2 , v, P\\ 2 , and P? h0 . Similarly, there are seven unknown 
variables in the cathode gas reaction region: C e , 4>, E c , Co 2 , t> , Po 2 » and The 

reactant gases in these layers will have the same flux expressions as developed in the 
gas diffusion regions. Since the gases dissolve into the electrolyte, the gas dissolution 
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Table I. Summary of governing equations for 
the anode gas diffusion region (z a < z < z a( i)- 
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Table II. Summary of governing equations for the 
cathode gas diffusion region ( z cd < z < z c ). 


d_( Pj dP °> \-n 

d£ y ^h 2 o 2(1 — fa)Po 2 J 


[ 32 ] 


P C T 


PO: 


+ P lhO 


d 2 E c 

d, e 2 


= 0 


[ 33 ] 


[ 34 ] 
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rates are included in the equation of continuity. Combining Eqs. [4], [5], and [13] 
gives an expression for the hydrogen gas pressure 




[35] 


in the anode and an expression for the oxygen gas pressure 

6j t = v • ~ ) 


[36] 


in the cathode. 

In order to account for the evaporation of water into the gaseous pores, linear 
correlations were developed to give the partial pressure of water as a function of 
KOH concentration and temperature (Appendix A): 


P\\ 2 0 = a H 2 0 + &n 2 0 • Ce 


[37] 


where an 2 o and ^n 2 o are regression coefficients at a constant temperature. Since the 
electrolyte concentration varies in the liquid pores of the gas reaction regions, the 
water vapor pressure will also vary as given by Eq. [37] and which will influence 
the gas phase transport as given by Eqs. [35] and [36]. Note that these correlations 
assume that equilibrium will be established between the gas and liquid phases in the 
gas reaction regions. 

The dissolved reactant gas concentrations are given by combining Eqs. [4], [6], 
[9], and [13] to give 


de'Cn, 

dt 


= V r ,V 2 Ch 2 - V 


(*’ Ch 2 ) 


lad 


a Ll a 

2 F 


+ a9 ap\ 1 2 


// h 2 ^ h 2 — C \\ 7 

8a 


[38] 


for the dissolved hydrogen concentration in the anode gas reaction region and 


dc l Ci 

dt 


o 2 


V' o y-Co, - V • (v Co,) +~^+ a'P'o, 


Ho, Po 2 — Co 2 


[39] 


for the dissolved oxygen concentration in the cathode gas reaction region. The local 
current densities in Eqs. [38] and [39] are expressed using the Butler- Volmer kinetic 
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expression, Eq. [11]: 


la = l. 


Ch.\ PH VG 


C H, 


c; 


exp 


>7 F 

“ “ {Ea-4>- Ua) 


— l 


•o 


exp - 


a a c n a F 

RT 


RT 

( E a -4>-U a ) 


[40] 


for the anodic current density and 

• - o \fCe 

lC C [Ue r 


a x? ' e 


exp 


C o 2 

c o 2 




exp 


<e£(e c - <p - u c ) 

E c -<t>- U c ) 


RT 


[41] 


for the cathodic current density. The superscribed r variables pertain to the reference 
concentrations associated with either the cathodic, i°, or anodic, i°, exchange current 
densities. The reference potentials, U, for each electrochemical reaction are given as 
a function of temperature, electrolyte concentration, and the partial pressure of the 
reactant gas relative to a standard reversible hydrogen electrode as defined by (55) 

V = U> - ^ Y, 1" (-) - E In f 42 l 

nf 4^ V P ) n RE F ^ V P J 

where ifi is the standard potential at temperature T and f/j[ E is the reference potential, 
both relative to the standard hydrogen electrode which is defined to be zero for 
convenience. Hence, these reference potentials are 


U c = Ui - M!±I ,n ( c ' 


for the cathode and 


Ua = K 


4 F 


8.314 T 

2F 


Po? 


(V) 


[43] 


In (Pu 2 Cl) (V) 


[44] 


for the anode. The standard electrode potentials at temperature T are given by (56) 
U* = 0.4011 - (T - 298.15) 1.682 x 10 -3 (V) [45] 
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and 

Ul = -0.828 - (T - 298.15) 8.360 x 10“ 4 (V) [46] 


The concentrations of the ionic species, K + and OH", are expressed by combining 
the material balance, Eq. [4], the ionic flux, Eq. [6], and the electrochemical rate, Eq. 
[9]. Combining these equations for each ionic species gives 

= X>+V 2 C+ + u+FV • (C+V<£) - V • (v m C+) [47] 

for the K + ions and 

= P_V 2 C_ u_FV • (C_V<£) - V ■ (tf C_) [48] 

for the OH - ions in both the anode and cathode gas reaction regions. In Eq. [48], 
the ia 1 term is given by i a a‘ a for the anodic gas reaction region and by i c a‘ c for the 
cathodic gas reaction region. The electroneutrality condition 

Y, «Ci = 0 [«] 

i 

can be used to relate the potassium and hydroxide ion concentrations in the solution 
phase, giving 

C' e = C_ = C+ [50] 

Equation [50] can be combined with Eqs. [47] and [48] to eliminate the K + and OH 
ion concentrations, giving: 

= I> + V 2 C e + u+FV • (C e V<f>) -V-(v m C e ) [51] 

and 

^ = D_V 2 Ce -y- U-FV • (CeV^) - V • (l>" C e ) 
for the anode and cathode gas reaction regions. 


[ 52 ] 



An expression for the volume average velocity can be formulated by first multi- 
plying the equation of continuity, Eq. [4], for each species / present in the solution 
by the partial molar volume (assumed constant) of each species 


Vi^- = [- V • N» + R-+ R?} Vi (i = +, o, 0 2 , H 2 ) 
and summed over all species resulting in 

d[e l (V-C~ + V+C+ + V 0 C 0 + V\\ 2 C\\ 2 + Vb 2 Co 2 )] _ 

dt 

-V • (N_V_ + N+V+ + N 0 V 0 + N 0a Vb 2 + N,i 2 Vji 2 ) 
+RLV- + V+ + R e 0 V 0 + Rq 2 Vq 2 4" R\ h hh 


[53] 


[54] 


Recognizing that 


1 = VLC_ + V+C+ + V 0 C 0 + V\hCu 2 + V 02 C 0a 


[55] 


and assuming that the partial molar volumes for oxygen and hydrogen are much smaller 
than for the electrolyte species, then Eq. [54] can be simplified to: 

() J 

— = -V • (N_VL + N+V+ + No Vo) 

+RIV- + R%V+ + R C 0 V 0 

The volume average velocity, v , defined by (52) is given as: 

v^Y^CiViVi = Y, N iVi 

i i 

which can be inserted into Eq. [56] and combined with Eq. [9] to give 

^ = -V • v - (s.V. + s+V+ + SoV 0 )^ [58] 


[56] 


[57] 


Newman (57) makes the following assumption for a binary electrolyte 



26 


which can be combined with 


1 =t°_ + t° + 

and the partial molar volume of the electrolyte 

K = u+V+ + v-V- 


to give 


and 


t°+V e = V-V- 


[60] 


[61] 


[62] 


t°_V e = i/_V+ [63] 

Combining Eqs. [60], [62], and [63] with the number of electrons transferred 

— n — + S-Z- [64] 

gives an expression relating the partial molar volume of the electrolyte to its dissociated 


ions 


s_V e t°_V e n 
+ s+V+ = + — 


[65] 


i/_ z. 


Since the liquid phase porosity does not change with time for the type of electrochem- 
ical reactions occurring in the alkaline fuel cell, the time rate of change of the liquid 
phase porosity is zero. Thus, combining Eqs. [58] and [65] gives an expression for 
the volume average velocity. 

's-V e . t°_V e n 


Vt = - 


+ 


V SoVo , r 

i/_~_ / nr 


i a 


[ 66 ] 


Since the total current density obtained from the fuel cell is equal to the integral 
of the local current densities as generated by the reaction, then 


l'R 

= -Jia‘ 


dz 


[67] 
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so that the potential drop is given by 

V 2 E a = — 

Oa 


[ 68 ] 


for the anodic gas reaction layer and by 


V' 2 E r 


[69] 


for the cathodic gas reaction layer where the anodic and cathodic current density 
expressions are given by Eqs. [40] and [41], respectively. 

Using Eq. [28], these governing equations for the anode gas reaction region are 
shown in their steady-state form in Table III. Similarly, the steady-state governing 
equations for the cathode gas reaction region are shown in Table IV. 


Separator 

In the middle of the alkaline fuel cell there is a porous separator matrix that is 
assumed to contain only solid non-conducting material and liquid— filled pores. The 
separator allows the ionic and dissolved reactant species to diffuse between the anode 
and cathode regions. The five unknown variables in this region are C e , <j>, v , Co 2 , 
and Ch 2 - The same development previously presented for the ionic species in the 
gas reaction regions can be applied here. No electrochemical reactions occur in the 
separator so that Eqs. [5 1 ] and [52] simplify to 

= P+ V 2 C e + u+FV • (C e V<f>) - V • (tf C e ) [86] 

and 

= l >_ V 2 C e - u-FV • (CeW) - v • (v C e ) = 0 [87] 

at ' ' 

Unreacted dissolved gases from the gas reaction regions may diffuse into the separator. 
Combining Eqs. [4] and [6] and simplifying gives 



[ 88 ] 



Table III. Summary of governing equations for 
the anode gas reaction region ( z a ,i < z < z aT ). 
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gL d ( 1%0 + P^ dPg 3 \ _ a9J)l 

RTL 2 dZ\P % 20 + faPu, d( 


a u \\? 


R\\ 2 F\\ 2 - Ch 2 


= 0 [70] 


P\uo = «h 2 o + % 2 o C e 


®fe a # 2 Cjl 2 1 ^ / ■ 

L 2 d £ 2 L dt V 


c Hl ) - 


7fn 2 Ai 2 - Ch 2 


= 0 


[71] 

[72] 


V + d 2 c e u+F 


L 2 dt 2 L 2 


d 2 <f> dC e d(j) 

^'e o ao 1 


Cedv V dC e 


d ( 2 dt dt\ L dt L dt 


= 0 


[73] 


V-d 2 C e ia a \i U -F 


L 2 dt 2 F I 2 


1 dv 


where 


d 2 <p dc e d<j> 
e dt 2 dt dt. 

s-V e t°_V e n a 


Cp dv V dC'e 


L dt L dt 


L dt 


+ 


+ s 0 V 0 


V- 


a'l 


a l a 


n a F 


1 d 2 E a i a a l a 


L 2 dt 2 



<»»/ 

RT 


( Ea-4-Ua ) 


0 [74] 


[75] 


[76] 


[ 77 ] 



Table IV. Summary of governing equations for 
the cathode gas reaction region ( z cr < z < z C( [). 
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PLo + Po 2 


dPo 2 


<V l o 2 


RT L 2 di\P{ h0 + 2(1 -f a )Po 2 d£ 

Pu 2 0 = a H 2 0 + b lh O C, 


Hq 2 Pq 2 - C o 2 


= ( 


go, d-Co , 

L 2 d( 2 ' 4 F 


+ 


?cflc + 


Ho 2 Po 2 — Co 2 


Cq 2 9v v dCo 2 


V+d 2 C e u+F 


d 2 <f> dc e d<j> 

e d( 2 di di 


L di L di 
C e dv v dC e 


L di L di 


= 0 


P_ d~C e l c a 'c 


where 


u_F 


L 2 di 2 F IF 


1 dv 


C t 


L di 


d 2 <t> 

dC e d<t>' 

'di 2 

+ di di\ 

S-V e 

t°_V e n c 

V- 

+ 


C e Ov V dC e 


L di L di 


= 0 


r e n c , ,7 a c l 

h -S 0 V 0 ) p 

z~ ) n c F 


1 d 2 E c i c a l c 


L 2 di 2 a c 


C \P' 




-i? 





exp 


a ^ cF (F c - * -U c ) 


RT 


3 [78] 

[79] 

0 [80] 
[81] 

[82] 

[83] 

[84] 


[ 85 ] 
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for the dissolved hydrogen gas and 

= V^V 2 Co, - v ■ (v C 0l ) [ 89 ] 

for the dissolved oxygen gas. The volume average velocity can be simplified from Eqs. 
[66] by recognizing that the local current density term is zero in the separator giving 

Vv" = 0 [90] 

These five governing equations for the separator region are summarized in Table V 
in their steady-state forms. 

C. Boundary Conditions 

There are six boundaries in the alkaline fuel cell that need to be incorporated into 
the model. Referring to Fig. 2, these boundaries occur at the following interfaces: 
anode gas channel/anode gas diffusion layer ( z a ), anode gas diffusion layer/anode 
gas reaction layer {z ad ), anode gas reaction layer/separator (z^), separator/cathode gas 
reaction layer (z cr ), cathode gas reaction layer/cathode gas diffusion layer (z cd ), cathode 
gas diffusion layer/cathode gas channel (z c ). Dirichlet or Neumann-type boundary 
conditions can be used to describe the phenomena that occur at these interfaces. These 
conditions will now be presented at each interface. 

Anode Gas Channel/Anode Gas Diffusion Interface 

The boundaries at this interface follow the Dirichlet-type conditions where the 
gaseous hydrogen and water vapor pressures are ‘fixed’ values. Since the total pressure 
at this interface remains constant, the partial pressures of hydrogen and water will vary 
so that the flux of water in the anode, Eq. [19], will be consistent with the fraction, 
f a , and the total flux of water predicted by the cell current density: 

r / 

N H 2 0 - 2f 


[ 96 ] 
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Table V. Summary of governing equations for the separator region ( z ar < z < z C r )• 


V+d 2 Ce u+F 

+ T 2 
L S L 


C 2 


d 2 (j> dc e d<j> 

v-e t 


Ce 5C e 


<9£ 2 Ls di Ls di 


= 0 


V- d 2 C e U-F 

4 oe 4 


d 2 <^ dc e d<j> 

a>o *T 


C e dv 


OCe 


de dZ d(\ Ls dZ Ls dZ 


= 0 


dv 


= o 


T>h 2 d-C)i 2 C h 2 dv v dCu 2 

L 2 S dZ 2 Ls dZ Ls dZ 


[91] 

[92] 

[93] 

[94] 


F{ h d 2 Co 2 Co 2 <9t>" u 5 Co 2 

If dZ 2 Ls dZ Ls dZ 


[ 95 ] 
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Hence, the water vapor pressure at this interface can be calculated by combining Eqs. 
[5], [14], [19], and [96] and solving for the water vapor pressure at the interface. Since 
the fraction of water leaving through the anode, f a , and the predicted cell current 
density, /, influence how much water is produced, the water vapor pressure at this 
interface is given as a function of f a and /. For the total pressure to remain constant, 
the hydrogen pressure is then given by 

Pth = p r ~ Pfho 1 97 1 

Since the model predicts the current for a given potential load, the anode potential at 
this interface can be arbitrarily set to any value. Hence, the anode potential, E a , was 
set to zero at this interface. These boundary conditions are summarized in Table VI. 

Anode Gas Diffusion/Anode Gas Reaction Interface 

At this interface, boundary conditions are needed for seven unknown variables: 
Ph 2 . P£ a0 . Ea. Ch», G' e , <f>, and v. The flux for the hydrogen gas as given by Eqs. 
[5] and [14] is continuous at this interface giving 

VttVPth \ B = v a ,vPth |„ I 10 '] 

The partial pressure of water at this interface is given by the correlation developed 
earlier for the partial pressure of water above a KOH electrolyte, Eq. [37]. The faradaic 
current density is continuous at this boundary as given by Eq. [25] resulting in 

* D VE a \ D =a R VE a \ R [ 102 ] 

The fluxes for the dissolved hydrogen and the K + and OH - ions are zero since 
there exists a solution/solid phase interface at this boundary. Setting the ionic flux 
expression, Eq. [6], to zero for each of these species results in: 

0 = VC'e 


R 


[103] 
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Table VI. Summary of boundary conditions at the anode 
gas channel/gas diffusion layer interface (z = z a ). 



PLo = f(faJ) 
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0 = Y<£ [104] 

R 

and 

0 = VC Il2 [105] 

The volume average velocity at this interface is given by Eq. [66] where the local 
current density is zero since no electrochemical reactions take place at this interface. 
Also, since this boundary is a solid/solution phase interface, the volume average 
velocity has to be zero. Hence, 

0 = Vt> and v = 0 [106] 

R 

These boundary conditions are summarized in Table VII. 
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Table VII. Summary of boundary conditions at the 
anode gas diffusion/gas reaction interface (z = z a d )• 
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Anode Gas Reaction/Separator Interface 

Eight boundary conditions need to be specified at this interface: for P\\ 2 , ^h 2 0’ 
E a , Ch 2 , C e , (f>, i", and C 0y Since hydrogen gas does not enter the separator, the 
flux of hydrogen gas, Eq. [5], is set to zero, giving 


VPh 2 



[114] 


and the partial pressure of water is again given by Eq. [37]. The faradaic current is 
zero at this boundary condition since the separator is non-conductive which simplifies 
from Eq. [25] to give 


V£ 


ft 



[115] 


The fluxes for the electrolyte species and dissolved hydrogen are continuous at this 
interface. Equating Eq. [6] as applied in the reaction layer and the separator layer 
for each species gives 


and 



p'b.VCh, 

= I’kVCn, 

s 

T> + VC c 

—v + FCeV<p 

= -v + vc e 

—n+FC e V<j) 

R 

R 

s 

-P-VG’e 

+u-FC e V<f) 

= -D_VCe 

+u-FC e V<t> 


R 

R 

s 


[116] 


[117] 


[118] 


The volume average velocity at this interface is given by equating Eq. [66] as applied 
in the anode gas reaction layer to that applied in the separator region. Note that 
the local current density term is zero since no electrochemical reactions occur at this 
boundary. 


Vo 


= Vo 


[119] 


R 


S 


Unreacted dissolved oxygen in the cathode gas reaction layer can diffuse through the 
separator towards the anode gas reaction layer. Any dissolved oxygen in the anode gas 
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reaction layer would be quickly consumed by an electrochemical reaction and would 
not significantly influence the system at all. Hence, the dissolved oxygen concentration 
can be set to zero at the anode reaction layer/separator interface. These boundary 
conditions for the anode reaction/separator interface are summarized in Table VIII. 

Separator/Cathode Gas Reaction Interface 

Similar conditions are used at this interface as used at the anode gas reaction 
layer/separator interface. Eight boundary conditions are needed to describe these 
variables: Pq 2 , P{ h0 , E c , C 0i , C e , <j > , ?>", and C\\ 2 . Since the phenomenon at this 
boundary condition are similar to that at the anode reaction/separator interface, only 
a summary of the boundary conditions will be presented as shown in Table IX. 

Cathode Gas Reaction/Cathode Gas Diffusion Interface 

The boundary conditions at this interface are analogous to those at the anode gas 
diffusion/anode gas reaction interface. Thus, only a summary of the seven required 
boundary conditions (Po 2 , P\\ 2 q-> E c , Cq 2 , C e , <t>, and v ) are presented in Table X. 

Cathode Gas Diffusion/Cathode Gas Channel Interface 

The boundary conditions at this interface are similar to those at the anode gas 
channel/anode gas diffusion interface except that the cathode potential is set to the 
applied cell potential. The partial pressures of oxygen and water are, again, allowed 
to vary at this boundary in an analogous manner as shown at the anode channel/gas 
diffusion interface. The boundary conditions for this interface are summarized in 


Table XI. 



Table VIII. Summary of boundary conditions at the 
anode gas reaction/separator interface ( z = z ar ). 
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[ 120 ] 

[ 121 ] 

[ 122 ] 
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0 = Co 2 
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[127] 
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Table IX. Summary of boundary conditions at the 
cathode separator/gas reaction interface (;? = z cr ). 




[ 129 ] 

[ 130 ] 


Po, ocq, _ go, ac 0 , 

Ls &Z s Lr dt ft 


V+ dC e 

u + FC e d<j> 

Is d£ 

1'S K 

V _ dCe 

u-FC e d<j> 

Ls dt 

S + Is H 


v + dc e 

u + FC e d<j> 

Lr (K 

ft L R dt 

V- dC e 

u-FCe d<f> 

Lft dt 

ft Lft dt, 


1 dv _ 1 dv 

T s ~di 5 = Tft~dl ft 


[ 131 ] 

[ 132 ] 

[ 133 ] 

[ 134 ] 




40 


Table X. Summary of boundary conditions at the 
cathode gas reaction/gas diffusion interface (z = z cd ). 


K, OPo , _ Vo, OPo, 

Lr r Ld d 

[136] 

P}\ 2 0 = a lhO + hhO Ce 

[137] 

(7 R 0E a <?D 9E a 
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Table XI. Summary of boundary conditions at the 
cathode gas diffusion/gas channel interface (z = z c ). 


Po, = Pi - PfhO [1431 

pk,o = f(f.,i) i 1441 


Ec — ^cell 


[145] 
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D. Model Parameters 


The alkaline fuel cell model requires various parameters that are specific to the 
cell structure, electrochemical reactions, and operating conditions. Table XII shows 
the electrochemical parameters used for the anode and cathode reactions for the 
base case conditions. The number of electrons transferred, charge numbers, and 
the stoichiometric coefficients for the cathode and anode reactions were taken from 
reactions [1] and [2], respectively, with the reactions following the form of Eq. [10]. 
Since the electrolyte is KOH, the dissociation constants of the ionic species are 1. The 
transfer coefficients for the hydrogen oxidation reaction in Table XII were estimated 
based on reported Tafel slopes of about 0.12 for the hydrogen evolution reaction (58). 
For the hydrogen evolution reaction, the Tafel slope, b, is 


b = 


2.3 


RT 

api a F 


[146] 


giving a“n a of about 0.5. Assuming that 


a„n -f a c n = n 


[147] 


then the product of the number of electrons transferred and the anodic transfer coeffi- 
cient would be 1.5 as shown in Table XII. Note that this assumes the mechanism for 
hydrogen oxidation and hydrogen reduction are the same. For the oxygen reduction 
reaction, the transfer coefficients are based on experimentally measured Tafel slopes 
of about 0.045 (59) giving the transfer coefficients shown in Table XII. Since reliable 
transfer coefficient data are not readily available for various temperatures, concentra- 
tions, potentials, and electrocatalysts, these transfer coefficients are assumed constant. 

Table XII also shows four adjustable parameters: the anodic and cathodic exchange 
transfer currents (?’° • «[, and i° ■ a £) and the anodic and cathodic diffusional film 
areas (a 9 a /S a and a%/8 c ). These parameters, as will be shown later, have a significant 
influence on the cell performance. For the base case conditions, these parameters were 
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Table XII. Electrochemical parameters for the 
anode and cathode reactions (base case conditions). 



Anode 


Cathode 

Parameter 

Value 

Parameter 

Value 

n a 

2 

n c 

4 

sh 2 

1 

so? 

-1 

S- 

2 

s __ 

4 

So 

-2 

^0 

-1 

z 

-1 


-1 

z + 

1 


1 

U- 

1 

V- 

1 


1 


1 

P h 2 

2 

<lo? 

0.5 

POH- 

2.0 

VOW- 

2.0 

K n a 

1.5 

( *a n c 

2.5 

ot a c n a 

0.5 

a c c n c 

1.5 

*« ■ a ‘a 

0.60 A/cm 3 

i°c ■ a ‘c 

0.60 A/cm 3 

a 9 a/S a 

5 x 10 8 cm' 2 

a g dS c 

5 x 10 8 cm" 2 



selected to give performance curves similar to those obtained experimentally. As will 
be shown later, these parameters can be adjusted to fit the model to experimental 
polarization data. 

Table XIII shows the structural and electrode parameters common to both the 
anode and cathode for the base case conditions. These parameters were selected as 
being representative of actual alkaline fuel cells. The base case operating conditions 
are shown in Table XIV where £ ce ii corresponds to the set cell potential. 

In order to obtain realistic performance predictions for the alkaline fuel cell, the 
diffusion rates of the various species were determined as a function of concentration 
(or pressure) and temperature. This was necessary to accommodate the changing 
pressure in the electrodes and the varying electrolyte concentration across the system. 
Appendix A shows the correlations for the diffusivity and solubility parameters. The 
gas phase diffusivities for oxygen in water and hydrogen in water are determined 
as functions of temperature and pressure based on corresponding states principles. 
The diffusivities of dissolved oxygen and dissolved hydrogen in KOH are based on 
experimental measurements at different temperatures and concentrations. The ionic 
diffusivities for K + and OH~ ions were scaled from their dilute solution values to 
concentrated solutions for different temperatures. Setschenow salt effect parameters 
were used to determined the solubilities of hydrogen and oxygen gas in KOH as 
a function of temperature and concentration in the form of Henry’s law constants. 
Based on these correlations, the reference concentrations necessary in the Butler- 
Volmer electrochemical reactions were determined using the base case conditions as 
shown in Table XV. Also shown in Table XV are the partial molar volumes of the 
electrolyte and the solvent (water). These partial molar volumes, assumed constant 
for this work, were obtained from (53) 



Table XIII. Structural and electrode parameters 
for the anode and cathode base case conditions. 
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Parameter 

Value 

Diffusion layer thickness, Lq 

0.0250 cm 

Reaction layer thickness, Lr 

0.0050 cm 

Separator thickness, Ls 

0.0050 cm 

Diffusion layer conductivity, ap 

5.0 S/cm 2 

Reaction layer conductivity, or 

5.0 S/cm 2 

Diffusion layer gas porosity, c 9 D 

0.70 

Reaction layer gas porosity, e g R 

0.10 

Reaction layer liquid porosity, e l R 

0.60 

Separator liquid porosity, t s 

0.80 

Cathode and anode tortuosity, r 

1.2 

Separator tortuosity, r 

1.0 
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Table XIV. Base case operating conditions. 


Parameter 

Value 

Initial electrolyte concentration, C e 

7 N 

Temperature, T 

O 

O 

O 

00 

Inlet anode gas pressure, Pj 

4.1 atm 

Inlet cathode gas pressure, Pf 

4.1 atm 

Applied anode potential, E a 

0.0 V 

Applied cathode potential, E c 

Ecell 
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for the partial molar volume of the electrolyte and 


V'o = 


M 0 


n dp 

p ~ c ezu: 


[149] 


for the partial molar volume of the solvent. The density of the electrolyte, as given 
by Akerlof and Bender (60), were correlated with temperature and concentration so 
that the change in density with concentration, dp/dC e , could be determined. Using 
the base case operating conditions, the partial molar volumes for the electrolyte and 
the solvent were then calculated using Eqs. [148] and 1149], respectively, as shown 


in Table XV. The transference number for the OH" ion with respect to the solvent 


velocity, t°_, introduced into the volume average velocity expression, Eq. [66], is 
determined by 


t° = 


A_ 


A_ -f A + 

where it is assumed constant over the electrolyte concentration range. 


[150] 


E. Method of Solution 

The alkaline fuel cell model consists of 25 governing equations and 36 outer 
and internal boundary conditions. These equations are highly coupled and non-linear 
so that a numerical method is necessary to solve the resulting system of equations. 
The model equations were discretized by using second order accurate finite— difference 
approximations (61,62) in the governing equation regimes and by first order accurate 
finite differences expressions at the boundaries. This switch from second order to first 
order accurate finite difference expressions at boundaries was necessary to keep a stable 
solution of the system of equations. The resulting finite difference approximations have 
a banded matrix structure that can be solved by Newman’s BAND(J) algorithm (53). 

As the alkaline fuel cell approaches the limiting current density, the dissolved 
oxygen concentration becomes prohibitively too small (~10 9 mol/cm 3 ) which creates 
numerical difficulties with other terms that are of a much higher order of magnitude. 
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Table XV. Reference concentrations and partial molar volumes. 
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The dissolved oxygen concentration had to be logarithmically transformed in order 
to achieve limiting current densities with the model. A logarithmically transformed 
variable, Cq 2 » can be defined by 


C* 

c o 2 


In 



[151] 


so that the derivatives become: 


PC p 2 

Oz 


Cb 2 


ex P \ c o 2 )~qT 


[152] 


for the first derivative and 


d 2 a 


o 2 


dz 2 


Cb,**P{Cch) 


dz 2 


+ 


dz 


[153] 


for the second derivative. Equations [152] and [153] can thus be substituted into the 
model equations for the dissolved oxygen concentration. 

To prevent forcing the electrolyte out of the separator, the total pressure drop across 
the separator was forced to zero by adjusting the fraction f a . Since a potentiostatic 
approach is used to model the alkaline fuel cell, the cell potential is set and the current 
density is calculated by the model. Convergency is obtained when the predicted cell 
current density is equal in all regions of the fuel cell as expressed by Eq. [25] for the 
electrode potentials in the gas diffusion regions, Eq. [67] for the integral of the local 
current densities in the gas reaction regions, and by 


I 

~nF 



( *— Oa , II2 ) 


[154] 


for the gas diffusion regions and 


I = F^ziNi (.= +.-) 


[155] 


for the separator. 
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F. Results and Discussion 

The base case parameters shown in Tables XII and XV were used with the alkaline 
fuel cell model to calculate the eleven dependent variables as a function of cell potential 
and spatial domain. The cell potential was varied over the potential range of 1.1 to 
0.7 V representative of the activation to concentration polarization range, respectively. 
These variables were then used to predict the current density as a function of cell 
potential. 


Dependent Variable Profiles 

To avoid forcing electrolyte out of the separator, the pressure drop across the 
separator region was forced to zero by adjusting the water vapor fraction, / a , until 
the pressures at the separator/reaction interfaces are equal as shown in Figs. 3 
and 4. These pressure distributions are influenced by the dissolution rate of the 
reactant gases and the evaporation rate of water. White et al. (32) showed that 
the electrode properties (e.g., a 9 ,a l ,c 9 ,e l ) were dependent on the system pressure 
for a pore spectrum that follows a Gaussian distribution. Since the variation of the 
anode and cathode pressures are small, the assumption of constant electrode properties 
throughout their domains is reasonable. For the base case conditions used in evaluating 
the model, the system reaches its limiting current density at cell potentials lower than 
about 0.85 V. At these potentials, the total pressures along the separator/reaction 
interfaces become constant at an appreciable pressure which indicates that gas-phase 
molecular diffusion does not become a limiting factor in obtaining the limiting current 
densities. A similar conclusion was found by (39) who determined, experimentally, 
that the gas permeability was not as important as the wettability of the catalyst in 
limiting the current density. The inclusion of Knudsen diffusion in the gas phase 
may result in significant voltage losses as shown by (28), however, the importance of 
Knudsen diffusion is still somewhat questionable as shown by (11). 
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Cathode 
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Figure 3. Total pressure distribution in the cathode. 
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Figure 4. Total pressure distribution in the anode. 
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Figure 5. Variation of the oxygen pressure in the cathode. 
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Figure 6. Variation of the hydrogen pressure in the anode. 
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Figure 7. Distribution of the water vapor pressure in the cathode. 
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Figure 8. Distribution of the water vapor pressure in the anode. 
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The total pressure profiles can be separated into their respective partial pressure 
profiles as shown for oxygen (Fig. 5), hydrogen (Fig. 6), cathode water vapor (Fig. 
7), and the anode water vapor (Fig. 8). As expected, the reactant gas pressures 
decrease very slightly in the diffusion and reactant layers at high cell potentials. 
However, at lower cell potentials, more of the reactant gases are consumed by the 
electrochemical reactions as evidenced by the steeper pressure drops. The reactant 
gas pressure profiles again confirm that at limiting current conditions the oxygen and 
hydrogen pressures are not a limiting factor since the gas pressures remain constant 
along the separator/reaction interfaces. The water vapor pressure profiles. Figs. 7 and 
8, show that more water evaporates into the anode reaction layer than the cathode 
reaction layer. 

The variation of the electrolyte concentration throughout the reaction layers and 
the separator is shown in Fig. 9 over the potential range of 0.7 to 1.1 V. Since the 
electrolyte does not circulate outside the fuel cell, the model assumes that no loss or 
production of the initial charge of electrolyte will occur. That is, the total number 
of moles of KOH is assumed to remain constant in the fuel cell. Previous models 
(23,28) for porous gas diffusion electrodes have typically neglected this electrolyte 
concentration variation. This simplification is reasonable at low current densities as 
shown by the relatively constant electrolyte concentration at high cell potentials in 
Fig. 9. However, at low cell potentials, an approximate 1.2 M change results in the 
electrolyte concentration from the anode to the cathode. This large variation has a 
significant impact on the evaporation of water and on the physical properties of the 
species present. Based on the equilibrium expression for the water vapor, Eq. [37], 
more water vapor will be present in the anode at low cell potentials than at high 
potentials. The accumulation of water vapor in the anode could flood the electrode at 
low cell potentials degrading the performance of the system. In the cathode regions, 
the larger electrolyte concentration at low cell potentials causes less water vapor to 



58 


be present which could dry out the electrode. Figure 9 also shows a fairly constant 
profile in the 0.7 to 0.8 V potential range indicating that the electrolyte concentration 
does not become a limiting factor in the limiting current region. 

The solubility of oxygen and hydrogen gas in KOH, as given in terms of Henry’s 
law constants, Eq. [171], are noticeably dependent on the electrolyte concentration as 
shown in Figs. 10 and 11, respectively. The dissolved oxygen concentration steadily 
decreases in the cathode reaction region as the cell potential is lowered obtaining an 
exceedingly small value (on the order of 1.0 x 10 7 M) at cell potentials lower than 
0.85 V. This decrease in the dissolved oxygen concentration results from an increasing 
consumption of dissolved oxygen by the electrochemical reaction, Eq. [41], and by 
a decrease in the solubility with increasing electrolyte concentration. Similar results 
occur for the dissolved hydrogen concentration as shown in Fig. 1 1 where the dissolved 
hydrogen concentration attains a steady value of about 0.36 x 10-3 M at cell potentials 
lower than 0.85 V. Comparing Figs. 10 and 11 shows that the low concentration of 
dissolved oxygen at potentials lower than 0.85 V causes mass transfer limitations 
in the cathode. Therefore, in order to obtain high current densities, more dissolved 
oxygen is needed in the cathode. The model of a single electrode by (28) predicts that 
the dissolved oxygen diffusion contributes the least amount to the polarization losses. 
However, this result was obtained for a low current density of 200 mA/cm 2 . At this 
current density, the AFC model predicts an appreciable amount of dissolved oxygen 
present in the electrolyte so that diffusional resistances of dissolved oxygen are not 
as significant as at the limiting current. 

The anode and cathode both experience potential drops during operation as shown 
in Figs. 12 and 13. These profiles show that the potential drops are typically less 
than 10 mV for electrodes with high conductivities of 5 S/cm. Combining these solid 
electrode potentials with the solution phase and reference potentials through Eq. [12] 
gives the overpotential or driving force for each electrode. As shown in Fig. 14, the 
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Figure 10. Dissolved oxygen concentration in the separator and cathode. 
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Figure 11. Dissolved hydrogen concentration in the anode and separator. 
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cathode requires a much larger overpotential than the anode at the same cell potential. 
The large overpotential in the cathode results from the low concentration of dissolved 
oxygen in the electrolyte. Attaining even larger overpotentials than shown in Fig. 14 
for the cathode would be difficult since the dissolved oxygen concentration is extremely 
low. The transfer rates associated with these overpotentials are shown in Fig. 15 where 
the transfer rates increase in magnitude from the gas diffusion/reaction interfaces to the 
separator. The location and distribution of the reaction zone has been investigated by 
(29,39) where it is has been concluded that most of the electrochemical reaction occurs 
within a 0.01 cm distance in the catalyst layer as measured from the gas diffusion layer. 
However, in addition to obtaining the optimal reaction layer thicknesses, it is equally 
important to determine the distribution of the current throughout the reaction regions. 
As shown in Fig. 15, the transfer currents have a significant variation across the 0.005 
cm thick reaction layers. Since the transfer currents are large near the separator, more 
electrocatalysts could be distributed near the diffusion/reaction interfaces to increase 
the reaction rates at these points. Optimizing the catalyst distribution with respect 
to the current density, amount of catalyst material, and cost could yield an improved 
electrode performance. 


Parameter Effects on Polarization 

For comparison purposes, the models’s predicted polarization is shown with some 
experimental polarization data (63) for two sets of operating conditions in Fig. 16. 
The parameters shown in Table XII were used for the low temperature and pressure 
polarization in Fig. 16. Close agreement is obtained between the predicted and ex- 
perimental polarization results. Since the complete operating conditions and fuel cell 
specifications are not given by (63) for the experimental results, an accurate compari- 
son between the model and experimental results cannot be made. Parameter estimation 
could be used to fit the model to the experimental data if more experimental infor- 
mation were known. This would allow the model to accurately predict the fuel cell 




















67 


performance beyond the experimental domain. Note that the model predictions for the 
high pressure and temperature polarization in Fig. 16 only cover part of the experi- 
mental data range. This is due to some of the inaccuracies of the transport property 
correlations at high temperatures and pressures. Obtaining reliable experimental data 
for the various transport properties at high temperatures and pressures would improve 
the model’s ability to predict high performance results for a wider range of operating 
conditions. 

The model predictions in Fig. 16 are influenced by four parameters: the anodic 
and cathodic exchange transfer currents (i° ■ a l a and i° ■ a l c ) and the anodic and cathodic 
diffusional film areas (a 9 u /8 a and a a c /8 c ), where the parameters shown in Table XII 
were used for the low temperature and pressure comparison. These parameters can 
significantly influence the polarization as shown in Fig. 17 for different values 
of the cathodic exchange transfer rate. This parameter is shown to influence the 
activation polarization region of the system without influencing the slope of the 
ohmic polarization region or the limiting current density at all. This suggests that 
an appreciable increase in the fuel cell performance up to the limiting current density 
can be obtained by increasing the electrocatalytic activity. Since the model predicts the 
same limiting current density for increasing exchange transfer currents, the dissolved 
oxygen concentration has to become even smaller to offset the higher exchange current 
densities. Similar results as shown in Fig. 17 were obtained for different anodic 
exchange transfer currents. 

The effects of different diffusional film areas are shown in Fig. 18 for the cathode 
and in Fig. 19 for the anode where the base case conditions were used. Two benefits 
are achieved by increasing the cathodic diffusional film area parameter. First, the slope 
of the ohmic polarization is minimized allowing larger current densities to be obtained 
and, second, the limiting current density is increased. Similar inferences also apply 
to the anodic diffusional film area parameter. Increasing a?/<5 c is associated directly 
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with increasing the oxygen gas dissolution rate as given by Eq. [13]. Hence, in order 
to obtain high current densities, the number of gas-liquid sites should be increased 
allowing more gas to dissolve into the electrolyte. Note that this does not necessarily 
imply that more gaseous-filled pores should be present in the three-phase electrode. 
That is, increasing the gas phase porosity in the reaction layer, t 9 R , will increase the 
number of gaseous filled pores, but not necessarily the gas phase specific surface area, 
a 8 . Doubling the anode gas phase diffusional film area from the base case value of 
5 x 10 8 cm -2 does not show any change from the base case polarization (shown by 
the solid line in Fig. 19). This simply indicates that the oxygen dissolution rate is 
limiting the performance. However, when a 9 a /8 a is lowered to 1.0 x 10 cm , the 
resulting performance decreases indicating that the hydrogen dissolution process is 
rate controlling. These results indicate that a significant interaction can occur between 
the anode and cathode in controlling the polarization. The ability to investigate the 
interaction of the anode, separator and cathode is one advantage to using a complete 
model of the alkaline fuel cell rather than using single electrode models as previously 
done. 

The performance of the alkaline fuel cell is also dependent on the values of the 
transfer coefficients used in the Butler- Volmer kinetic expressions, Eqs. [40] and [41] 
as shown in Fig. 20 for different anodic and cathodic transfer coefficients. These 
parameters are shown to influence the polarization in an analogous manner as the 
cathodic and anodic exchange transfer currents. The determination of these transfer 
coefficients is important in order to obtain better model predictions. 

Effects of Operating Conditions 

The effect of increasing the total pressure in both the anode and cathode gas 
channels would increase the performance of the fuel cell as shown in Fig. 21 for 
three different pressures. Since earlier results showed that the dissolved oxygen 
concentration had a limiting effect on the performance, increasing the pressure allows 
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Figure 18. Influence of the cathodic diffusion film area on the 
fuel cell polarization evaluated under the base case conditions. 
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Figure 19. Influence of the anodic diffusion film area on the 
fuel cell polarization evaluated under the base case conditions. 
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more gas to dissolve into the electrolyte as given by Henry’s law, Eqs. [171] and 
[172], which increases the fuel cell performance. The model predictions give no 
apparent upper bound on the system pressure where the performance would start to 
degrade due to transport or kinetic effects. 

Figure 22 shows the results of different system temperatures on the alkaline fuel 
cell performance. Increasing the temperature is shown to increase the activation 
polarization while extending the concentration polarization limit. As the system 
temperature is decreased, the limiting current density increases while the cell potential 
at the limiting current decreases. These results suggest that an optimal temperature 
might exist that is higher than 373 K that would give a maximum in the power density. 
Since the temperature has a significant impact on the cell performance, the inclusion 
of a thermal balance into a future alkaline fuel cell model could improve the model 
predictions. 

The effects of different average electrolyte concentrations on the fuel cell perfor- 
mance are shown in Fig. 23 where a decreasing electrolyte concentration increases 
the performance. The average electrolyte concentration in most alkaline fuel cells is 
about 7 N (32%) since this concentration corresponds to the highest conductivity of the 
electrolyte. The model results in Fig. 23 apparently contradict this fact. The model 
presented here indicates that the dissolved oxygen concentration is a major factor in 
determining the fuel cell performance, based on the operating and cell parameters 
shown in Tables XII-XIV. Since the solubility of oxygen increases with a lowering of 
the electrolyte concentration, Eq. [171], a higher cell performance results when more 
oxygen dissolves into the electrolyte. 

In order to determine the major limitations to the alkaline fuel cell performance, 
the base case conditions were used to examine the performance when certain forms 
of resistance were neglected as shown in Fig. 24 where the percent increase in 
the predicted current densities over the base case current densities are shown. To 
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Figure 24. Improvement in the current density with respect to the base case 
current density x 100) when different forms of resistance are neglected. 
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accomplish this, the various transport parameters were set to large values such that a 
further increase in the transport parameter would not change the resulting polarization. 
Thus, to investigate the effects of no gas phase resistances, the gaseous diffusivities 
were set to large values (~ 10000) so that no gas phase transport would exist. In 
an analogous manner, other transport parameters were increased to minimize their 
respective form of resistance. As shown in Fig. 24, minimizing gas phase diffusional 
resistances will contribute the least to improving the performance. Ionic resistance 
effects were minimized by increasing the conductivity of the electrolyte which gave 
a better performance increase than neglecting gas phase diffusional resistances. The 
common assumption of no electronic drop is shown in Fig. 24 to give the highest 
increase in performance for cell potentials greater than 0.9 V. However, for cell 
potentials lower than 0.9 V, the solution phase resistances have the most influence 
on the performance. This is reasonable since at potentials that approach the limiting 
current density, mass transfer effects prevent the attainment of even higher current 
densities. Increasing the liquid phase diffusion rates and the solubilities of the reactant 
gases will increase the alkaline fuel cell performance at low cell potentials. 
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CHAPTER IV 

PARAMETER SENSITIVITY AND OPTIMIZATION 
PREDICTIONS OF THE ALKALINE FUEL CELL 

A. Introduction 

In designing high performance alkaline fuel cells, there are various attributes 
that can significantly influence the system. Such attributes might be the gas and 
liquid phase porosities, reaction layer and separator thicknesses, or the number of 
gas-liquid sites in the three phase electrodes. One way to investigate the relative 
importance of these parameters is to use a mathematical model that describes the 
chemical, electrochemical, and physical processes occurring in the fuel cell. Typically, 
models of single, three phase electrodes are used in determining polarization losses 
and optimal design parameters (28,64). However, these models do not consider 
any interactions between the anode, cathode, and separator which can significantly 
alter the performance of the system as well as alter the optimal values for certain 
parameters. A sensitivity analysis is performed on various parameters to determine 
which parameters are the most influential in increasing or decreasing the current 
density. This information can indicate the direction one should take in order to 
design better fuel cells. The results of the sensitivity analysis can also suggest which 
parameters should be obtained with more accuracy through further modeling studies 
or through experimentation. 

To achieve high performance in the alkaline fuel cell, various design parameters 
can be optimized so that the fuel cell will deliver the maximum attainable power 
density. The important design parameters in the alkaline fuel cell are the thicknesses 
of the anode and cathode diffusion and reaction layers (Lp, Lr), separator thickness 
( Ls ), electrode conductivity (a), gas and liquid phase porosities ( e 8 , e l ), and the gas 
and liquid phase specific surface areas (at, d). By using the detailed model of the 
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alkaline fuel cell, these parameters are investigated in order to determine if an optimal 
value exists for each parameter. The sensitivity of the model predictions to various 
parameters will be examined next followed by the determination of the optimal design 
parameters to maximize the alkaline fuel cell’s power density. 


B. Sensitivity Analysis 

In order to determine the relative importance of the transport, kinetic, and structural 
parameters on the fuel cell’s performance, a sensitivity analysis can be applied. The 
sensitivity analysis can indicate which parameters have the largest influence on the 
predicted current density and, also, over which range of cell potentials the parameters 
have the most influence. Additionally, the sensitivity analysis can indicate which 
parameters are capable of being estimated when the model is used in conjunction 
with experimental data and a parameter estimation technique. That is, if a small 
perturbation in a parameter does not significantly change the predicted current density, 
then that parameter could assume a large range of values, all of which will give the 
same performance. The sensitivity coefficient can be defined as the difference in the 
predicted current density to a base-case current density for a small, dimensionless 
perturbation in a parameter j, while holding all other parameters constant: 


dl 

d in Oj 


i - r 


(A /cm 2 ) 


[ 156 ] 


where 0* and I* are the base-case parameter value and current density, respectively. 
Hence, large sensitivity coefficients indicate the parameter of interest significantly 
influences the current density. Large sensitivity coefficients may also indicate which 
parameters should be obtained with more accuracy through further modeling or 
experimental studies. That is, if the value for a parameter is not accurately known 
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and the parameter has a large sensitivity coefficient, then that parameter value should 
be ascertained with more accuracy to gain confidence in the model predictions. 

All sensitivity coefficients calculated for this work were accomplished by increas- 
ing the parameter of interest by 5 % over the base case values (shown in Tables XII — 
XV) and calculating the resulting change in the current density as given by Eq. [156]. 
This was performed over the potential ranges of (0.8 — 0.85 V), (0.85 — 0.93 V), and 
(0.93 - 1.1 V) representative of the concentration, ohmic, and activation polarization 
regions, respectively, for the conditions of the fuel cell simulation. By calculating 
the sensitivity coefficients in this manner, the effects of the parameters on the current 
density can be investigated under conditions of the various forms of polarization. 

Sensitivity of Transport and Electrokinetic Parameters 

The sensitivity coefficients for various parameters are shown in Figs. 25 and 
26 for parameters specific to the cathode and anode, respectively. As shown in 
the activation and ohmic regions, the model predictions are most sensitive to the 
transfer coefficients, liquid phase specific surface area, and the reactant gas reaction 
orders. Since the model predictions are extremely sensitive to the transfer coefficients 
as governed by the exponential terms in the Butler- Volmer expression, Eq. [11]» 
small perturbations in the transfer coefficients can significantly affect the predicted 
current density. Unfortunately, obtaining accurate values for the transfer coefficients 
is difficult since they vary too much depending on the temperature, cell potential, 
electrocatalyst, and electrode structure. Parameter estimation techniques could be used 
to fit the model to reliable experimental data by predicting the values for the transfer 
coefficients. This may necessitate a reformulation of the Butler- Volmer expression as 
shown by (65) to avoid numerical difficulties in the parameter estimation method. 

The model predictions show little sensitivity to small perturbations in the con- 
ductivities of the cathode and anode diffusion regions as shown in Figs. 25 and 26, 
respectively. The dependence of the model predictions on the conductivities in the 
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Figure 25. Sensitivity of the model predictions to cathode parameters for 
the concentration (I), ohmic (II), and activation (III) polarization regions. 





Figure 26. Sensitivity of the model predictions to anode parameters for 


the concentration (I), ohmic (II), and activation (III) polarization regions. 
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reaction layers was even less pronounced than in the diffusion layers. The relative 
insensitivity of a 9 a / 8 a on the model predictions over the entire range of cell potentials 
indicate that the dissolution of hydrogen gas into the electrolyte is not rate limiting. 
The most influential parameter in the concentration polarization region is a 9 /8 c which 
governs how much oxygen gas dissolves into the electrolyte through Eq. [13]. 

Sensitivity of Thickness Parameters 

The effects of fuel cell thickness on the predicted performance are shown in Fig. 
27 where the cathode reaction layer thickness is the most sensitive to the model 
predictions, followed by the anode reaction layer and separator thickness. The anode 
and cathode gas diffusion layer thicknesses are shown to have little effect on the 
model predictions. It is also apparent in Fig. 27 that the limiting current density can 
be increased by increasing the cathode reaction layer thickness or by decreasing the 
separator thickness. However, as will be shown later, increasing the cathode reaction 
layer thickness too much can degrade the performance. 

Sensitivity of Porosity Parameters 

The effects of porosity on the model predictions are shown in Fig. 28 where 
e l R c has the largest influence on the model predictions followed by e 9 Rc , e s , and t l R a . 
Increasing the gas phase porosities in the diffusion layers and in the anode reaction 
layer showed no change in the model predictions. Since the diffusivities are influenced 
by the porosities as indicated by Eq. [7], Fig. 28 suggests that diffusion is not as 
important at high cell potentials as at low cell potentials. This is verified by Fig. 29 
showing that the diffusion coefficients have the most influence on the cell performance 
in the ohmic and concentration polarization regions. 

C. Current Density Optimization 

Thickness Effects 

The previous analysis on the sensitivity coefficients showed that small perturba- 
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tions in design parameters could yield significant improvements in the current density. 
However, the sensitivity analysis does not allow a quantitative prediction on what 
values the design parameters should have in order to provide the best performance. 
By using the mathematical model of the alkaline fuel cell, various design parameters 
can be optimized so that the system achieves the maximum attainable power density. 
From Fig. 27, it can be seen that the anode and cathode reaction layer thicknesses 
and the separator thickness have the most effect on the performance. Calculating the 
limiting current density for these parameters as they are varied individually over a 
20 to 300 micron range with the others set equal to their base case values (Tables 
XII - XV) gives the results shown in Fig. 30. As shown in Fig. 30, a maximum 
occurs in the limiting current density for L a R at about 40 fi m and for L C R at about 
225 fj,m. The separator thickness does not show a maximum in the limiting current 
density indicating that its thickness should be as small as possible. Kenjo and Kawatsu 
(39) measured a flat limiting current density of about 1.5 A/cm 2 corresponding to a 
reaction layer thickness of 100 to 270 /xm for an oxygen electrode. Although different 
operating conditions were used, the location of the optimal thickness range in Fig. 30 
for L c r is similar to that obtained by (39). It has commonly been thought that in- 
creasing the reaction layer thickness should increase the limiting current density since 
more reaction sites are present in the electrode. However, according to our model, 
the reason for the decrease in the limiting current density beyond an optimal thickness 
is due to a lowering of the gas solubility, not to gas phase diffusional resistances. 
Increasing the reaction layer thickness causes a higher electrolyte concentration in the 
cathode which lowers the solubility of the reactant gas into the electrolyte. This will 
diminish the liquid phase diffusion of the dissolved reactant gas leading to lower fuel 
cell performance. 

Porosity Effects 

The sensitivity analysis for the porosity parameters, Fig. 28, show that the 

0 .- 3 . 
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concentration and ohmic polarization regions are significantly influenced by the various 
porosities. To investigate the optimal values for these porosities, the limiting current 
density was calculated for different parameter settings for the liquid phase porosities. 
These results are shown in Fig. 31. Note that the anode and cathode reaction layers 
were assumed to have a total porosity of 0.7 causing a constraint for the gas and liquid 
phase porosities 


0-7 = <4 + <4 


[157] 


As shown in Fig. 31, increasing e l R a does not cause any noticeable difference in the 
predicted limiting current density. However, in the cathode, a dramatic increase in the 
limiting current density occurs up to an optimal porosity of about 0.695 where upon 
further increasing the porosity a rapid decline in the limiting current density occurs. 
Since the model treats the gas phase transport as occurring by molecular diffusion and 
convection, only a small fraction of the total electrode porosity is needed for the gas 
phase. As the gas phase porosity approaches a small value ( e.g ., 0.005), the effective 
gas phase diffusivity becomes even smaller through Eq. [7] resulting in mass transfer 
limitations for the gas phase transport. Increasing the separator porosity increases 
the limiting current density more rapidly at lower porosities than at higher porosities. 
Since the separator porosity was assumed to be 0.8, further increasing the porosity will 
result in a slight increase of 30 mA/cm 2 at the limiting current density. Unfortunately, 
a maximum does not occur in the limiting current density for the separator porosity 
preventing an optimal porosity from being recognized. Note that other criteria could 
be considered in determining an optimal separator porosity in addition to an extremum 
in the current density such as the mechanical strength of the separator or lifetime, but 
this is beyond the scope of this work. 



(gino/v) 
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Specific Surface Area Effects 

The effects of the three specific surface areas ( a? a , a?, a l c ) on the limiting current 
density are shown in Fig. 32. The performance curve for a} a is similar to that 
obtained for a} c and, thus, is not shown in Fig. 32. Note that since the diffusion 
layer thickness, 8, and the exchange current density, i° , always appear with either a 
gas or liquid phase specific surface area, 8 was arbitrarily set to 1.0 x 10 -5 cm and 
i° to 1.0 x 10 -5 A/cm 2 . It should also be noted that a l a and a[ do not affect the 
limiting current density as shown in Fig. 17, but they do affect the polarization prior 
to reaching the limiting current density. Thus, to properly investigate the effects of 
these parameters, the current density at 0.9 V was predicted for a l c ranging from 500 
to 50,000 cm 2 /cm 3 as shown in Fig. 32. The predicted current density approaches 
a fairly constant value as a[ is increased, suggesting that further increases in a* will 
only contribute a marginal improvement to the current density or that a limiting current 
density has been attained. In Fig. 32, it is shown that increasing the anode gas phase 
specific surface area beyond 2000 cm 2 /cm 3 does not improve the current density at 
all. Although increasing this parameter increases the dissolution rate of hydrogen gas 
into the electrolyte as governed by Eq. [13], the dissolution rate for oxygen is still 
limiting the current density. This is verified by Fig. 32 for a 9 where increasing this 
parameter causes an increase in the current density. The large increase in performance 
due to increasing a 9 shows the importance of designing three-phase electrodes with 
large interfacial gas-liquid surface areas, especially for the oxygen electrode. Since 
the dissolution rate is dependent on the gas phase specific surface area, an increase in 
a 9 c will allow more gas to dissolve into the electrolyte to react. 

Optimal Power Density 

To achieve optimal performance in the fuel cell, the more influential parameters 
of the model can be used to maximize the predicted current density or power density. 
It should be noted that other criteria than a maximum current density could be used 




a x 10 -3 (cm 2 /cm 3 ) 

Figure 32. Effects of the gas and liquid phase 

specific surface areas on the current density. 
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in formulating an objective function. For example, Newman (66) optimizes an acid 
fuel cell by considering the average current density and the utilization of hydrogen 
based on capital, power, and fuel costs. Since the main objectives of this work are to 
increase the maximum attainable power density, economic factors are not considered. 
To maximize the power density, an objective function can be defined as: 

max P0 = £cell ■ I (&, £cell) [158] 

where 0 represents a vector of unknown parameters 9j and £ceii is itself an unknown 
cell potential. The optimal parameters and cell potential can be selected such that the 
power density as given by Eq. [158] is at a maximum. 

The results of the sensitivity analysis shown in Figs. 25 - 28 and the single 
parameter optimal studies shown in Figs. 30 - 32 indicate which parameters can be 
optimized in Eq. [158]. The anode and cathode reaction layer thicknesses and the 
liquid phase porosity in the cathode reaction layer are the only parameters that caused 
an extremum in the predicted current density. All other parameters investigated here 
monotonically increased the current density as the parameters where lowered ( e.g.,Ls ) 
or increased (<?.g., a®, ( s ) • Note that the transfer coefficients could be included 

in the optimization procedure since they have a strong effect on the current density 
as shown in Figs. 25 and 26. Additionally, the reactant gas reaction orders, qo 2 and 
pn 2 could also be included in the optimization procedure. However, since it may be 
more difficult to manufacture electrodes with certain values for the transfer coefficients 
and reaction orders than for the thicknesses or porosities, the transfer coefficients and 
reaction orders were not considered in the optimization procedure. 

The IMSL routine UMINF (67) was used to minimize Eq. [158] by using a 
quasi-Newton method to determine L a R , L C R , e l R c , and £ceii- The optimized values are 
shown in Table XVI along with their starting values. The optimal cathode reaction 
layer thickness in Table XVI corresponds to about the same optimal value as shown 
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in Fig. 30, whereas a large difference results for the optimal anode reaction layer 
thickness. Since increasing L R above 55 j/m did not improve the current density as 
shown in Fig. 30, L a R was scaled over a thickness range of 50 to 300 //m using 
the optimal values in Table XVI for L C R , e l Rc , E ce n to investigate whether L a R is 
indeed at an optimal value. The resulting performance curve verified that the optimal 
anode reaction layer thickness shown in Table XVI does cause an extremum in the 
power density. Using the optimal parameter values other than E ce ii, the performance 
of the fuel cell was predicted by the model as shown in Fig. 33 in comparison to the 
base case performance. As can be seen, a significant improvement in the maximum 
attainable power density has been achieved just by optimizing the cell potential and 
three design parameters. 
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Table XVI. Optimal parameter values for maximizing the power density. 


Starting Values 
L\ = 0.005 cm 
L c r = 0.005 cm 

4,c - 0-«5 

£«ll = 0.75 V 


Optimized Values 
L\ = 0.01627 cm 
L% = 0.02234 cm 

4.c - o - 674 

£„ll = 0.803 V 
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CHAPTER V 

CONCLUSIONS AND RECOMMENDATIONS 

A. Alkaline Fuel Cell Model Conclusions 

A mathematical model of an alkaline fuel cell has been developed to predict the 
performance of the system for different cell potentials operating under steady-state 
and isothermal conditions. The model describes the phenomena occurring in the gas, 
liquid, and solid phases of the anode and cathode gas diffusion regions, the anode 
and cathode reaction layers, and the separator. The model accounts for the one- 
dimensional transport of reactant gases, water vapor, solution phase concentrations, 
solid and solution potential variations, and the volume average velocity. Performance 
results were obtained for a set of base-case conditions that could be used for a high 
performance alkaline fuel cell. Gas phase diffusional resistances were found not to 
significantly influence the performance of the system. The model predicts that the 
diffusion of dissolved oxygen contributes the most to the polarization losses at low 
potentials while the electronic resistances contribute the most resistance at high cell 
potentials. To obtain better performance with the alkaline fuel cell, it is suggested that 
the three-phase electrodes should be highly conductive and contain a large number of 
gas-liquid sites allowing more reactant gas to dissolve into the electrolyte. Increasing 
the pressure and temperature will also result in improved performance. 

Various attributes of the system such as the catalyst distribution, exchange transfer 
currents, and diffusional film areas could be optimized to yield better performance. 
The model shows that interactions between the anode and cathode exist suggesting 
that models or experiments based on full-cells are necessary instead of half-cells 
when estimating unknown parameters or optimizing various attributes. One of the 
main advantages of this model is its ability to quantitatively show the influence of 
different parameters on the predicted current density. These quantitative results can 
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help in the design of alkaline fuel cells as well as help focus the direction of future 
research on alkaline fuel cells. 

A sensitivity analysis of an alkaline fuel cell model indicates that many parameters 
can significantly influence the performance of the system, especially in the ohmic 
and concentration polarization regions. In particular, parameters specific to the 
oxygen electrode such as the reaction layer thickness, liquid phase porosity, gas 
phase specific surface area, and the cathodic transfer coefficient have been found 
to influence significantly the performance. The effect of various design parameters 
on the limiting current density have been investigated to determine if optimal values 
exist for the parameters. The model has shown that the anode and cathode reaction 
layer thicknesses, the liquid phase porosity in the cathode reaction region, and the 
cell potential can be optimized to give the maximum attainable power density. The 
optimal reaction layer thicknesses are shown to be a compromise between the number 
of reaction sites and the solubility of the reactant gases. A small fraction of the total 
porosity in the cathode reaction region is needed in the gas phase to sustain a high gas 
phase diffusion rate while maintaining a large diffusion rate in the liquid phase. The 
model predictions indicate that the largest improvement in the fuel cell performance 
will be recognized by increasing the gas phase specific surface area in the cathode 
followed by increasing the electrocatalytic activity or liquid phase specific surface area 
and decreasing the separator thickness from the base case conditions. 

B. Recommendations for Further Studies 

Future studies on the alkaline fuel cell should include a closer comparison of 
the model and experimental data. Incorporating more features ( e.g ., teflon content, 
catalyst distribution, pore sizes, etc.) of experimental three phase electrodes into the 
model may give a better description of the fuel cell performance. It may also be 
desirable to include air as the reactant gas in the oxygen electrode rather than oxygen 
itself since some applications of the alkaline fuel cell use air as the reactant gas. This 
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would require the addition of another reaction in the system since the C0 2 in the air 
would react with the OH~ ions causing the COjf to precipitate in the pores and the 
OH - concentration to be reduced. 

One desired extension to the alkaline fuel cell model would be to include the gas 
channel effects. As the reactant gases flow through the channels, their pressures 
decrease while the water vapor pressures increase in the anode and cathode gas 
channels. If too much water is produced, the electrodes may become flooded causing 
severe performance degradation. Including the gas flows in the gas channels may help 
determine when conditions of flooding will occur and suggest methods to minimize 
the effects of flooding. 

Another important effect to consider is the thermal management of the fuel cell. 
The heat produced by the alkaline fuel cell reactions is typically removed by a 
circulating electrolyte. However, when a non-circulating electrolyte is used, as in 
this work, severe thermal gradients may occur in the electrodes and separator which 
can further degrade the performance of the fuel cell. Including a thermal balance 
would allow an even more realistic description of the fuel cell. 



NOMENCLATURE 


Roman Symbols 

a g Specific gas phase surface area, cm 2 /cm 3 

a 1 Specific liquid phase surface area, cm 2 /cm 3 

b Tafel slope, V/decade 

C, Concentration of species t, mol/cm 3 


Di 

Vi 

E 

fa 

fc 

F 

H 


Hi 





n 


Free-stream diffusivity of species t, cm 2 /s 
Effective diffusivity of species i, cm 2 /s 
Electrode potential, V 

Fraction of water generated that leaves through the anode 
Fraction of water generated that leaves through the cathode 
Faraday’s constant, 96,487 C/mol 
Hessian matrix 

Henry’s constant for species i, mol/(cm 3 atm) 

Local current density. A/cm 2 
Exchange current density, A/cm 2 
Total cell current density, A/cm 2 
Limiting current density. A/cm 2 
Setschenow salt effect parameter, cm 3 /mol 
Thickness, cm 

Molecular weight of species i 
Number of electrons transferred 


Ni 


Pi 

P 


Pi 


Qi 

R 


Flux of species i, mol/(cm 2 -s) 

Anodic reaction order for species i 

Power density, W/cm 2 

Pressure of species i, atm 

Cathodic reaction order for species i 

Gas constant, 8.314 J/(mol-K) or 82.057 cm 3 -atm/(mol-K) 
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R* 

RP i 

Si 

t 

*? 

T 


U 

u e 


Ui 


v 

Vi 


Xi 

yi 

z 

Zi 


Electrochemical reaction rate, mol/(cm 3 -s) 

Transport rate across phase boundary, mol/(cm 3 -s) 

Stoichiometric coefficient of species i 
Time, s 

Transferrence number of species i relative to the solvent velocity 
Temperature, K 

Theoretical open-circuit potential evaluated at reference 

concentrations, V 

Standard electrode potential, V 

Reference electrode potential, V 

Mobility of species /, mol-cm 2 /(J-s) 

Volume average velocity, cm/s 

Partial molar volume of species i, cm 3 /mol 

Liquid phase mole fraction of species i 

Liquid phase mole fraction of species i in pure water 

Vapor phase mole fraction of species i 

Spatial coordinate, cm 

Charge number of species i 


Greek Symbols 


a a 

OL c 

8 

e 

V 

e 

°j 

A 


A, 


Anodic transfer coefficient 
Cathodic transfer coefficient 
Diffusion layer thickness, cm 
Porosity 

Overpotential, V 
Vector of parameters, 0j 

Parameter j 

Equivalent conductance, cm 2 /fl 

Limiting ionic conductivity of species i, cm 2 /ft 

Ionic conductivity of species i, cm 2 /0 
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Number of cations or anions produced by the dissociating 
electrolyte 

Dimensionless spatial coordinate 
Electrolyte density, g/cm 3 
Electrode conductivity. S/cm 2 
Tortuosity 

Solution phase potential, V 


Superscripts and Subscripts 

a Anode 


Cathode 
Diffusion layer 

Electrolyte or electrochemical rate 

Gas phase 

Species i 

Species j 

Liquid phase 

Reference condition 

Reaction layer 

Production rate due to transport across a phase boundary 
Separator layer 
Total value 
Cation (K + ) 

Anion (OH ) 

Solvent (water) 

Base-case value or logarithmically transformed variable 
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APPENDIX A 

CORRELATIONS FOR DIFFUSIVITY 
AND SOLUBILITY PARAMETERS 


In order to obtain a realistic simulation of the alkaline fuel cell, the various 
transport and thermodynamic properties should be dependent on the temperature 
and concentration since the concentration of all species readily changes during the 
simulation. Correlations were developed to predict the values of these transport and 
thermodynamic parameters as a function of concentration (or pressure) and the system 
temperature. 

Correlation of the Water Vapor Pressure to the Electrolyte Concentration 

Since it is assumed that equilibrium is established at the gas-liquid interfaces in 
the active regions of the fuel cell, the water vapor pressure above the KOH electrolyte 
can be correlated with the KOH concentration. Experimental measurements of this 
water vapor pressure above KOH have been measured from 20 C to 80 C for KOH 
concentrations ranging from 0% to 50% KOH (by weight) (76). Fitting this data to a 
linear function (R 2 » .99) allows the following correlation to be defined: 

Pe 2 o = «h 2 o + &h 2 o C t t 159 l 

where C e is in mol/cm 3 and the vapor pressure of water is in atm. The regression 
coefficients are shown in Table XVII. Note that for temperatures higher than 80 °C 
the water vapor pressure was extrapolated from the lower temperature correlations. 
Extrapolating the data in this manner can introduce some uncertainty in the water 
vapor pressure and, thus, in the model predictions. Hence, caution should be used for 
temperatures higher than 80 °C. 



Table XVII. Regression coefficients for the correlation of 
the water vapor pressure to the electrolyte concentration. 


Temperature (K) 

m 2 o 

& h 2 o 

333 

0.19503 

-11.222 

343 

0.30777 

-18.468 

353 

0.46858 

-28.095 
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Correlations for Diffusivities 


The gas phase diffusivities for oxygen in water and hydrogen in water, at low 
pressures, can be given by corresponding states principles (52) where 


2.334 


PD 


•» 


(Pc.P c y / HTc<T' l f n Jriir l + uw; 


= 3.640 x 10 


-4 1 


y/TJZ' 


[160] 


Equation [160] simplifies to 



4.2076 x 10 -7 r 2 334 

P 



[161] 


for oxygen in water and 



2.1410 x 10- 6 T 2 - 334 
P 



[162] 


for hydrogen in water where T is in °K, and P in atm. 

The diffusivity for dissolved oxygen in KOH and for dissolved hydrogen in KOH 
was obtained from the experimental measurements of Tham et al. (77) where they 
measured the diffusivities from 25 °C to 100 °C for 0 to 50% KOH solutions. These 
diffusivities were correlated as a function of KOH concentration by fitting their data 
to cubic polynomials. For temperatures higher than 80 °C, the data were extrapolated 
from the lower temperatures. 

The diffusivity of the electrolyte species were determined by scaling the equivalent 
conductances of the KOH electrolyte at 25 °C to the temperature of interest by using 
the ratio of the limiting equivalent ionic conductivities. That is, the limiting ionic 
conductivities given by (78) for 25 °C and 100 °C, shown in Table XVIII, were 
interpolated to the temperature of interest as shown in Table XVIII, for example, at 
80 °C. The ratio of the limiting ionic equivalent conductivity for species i, A°, to the 
limiting equivalent conductivity, A 0 (= A°_ + A^) , at the temperature of interest was 
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then multiplied by the equivalent conductivity, A, of the KOH electrolyte (at 25 C) 
as given by (79) to get the ionic equivalent conductivity. A*, for species i. 

A = 272.0 - 3846.59 <7 e 05 + 1.41 x 10 5 C e ( 1.0 - 7.191 C e 0 ' 5 ) (cm 2 /0) [163] 


where the electrolyte concentration is in mol/cm 3 . For example, at 80 C and for a 
7M KOH electrolyte, the ionic equivalent conductivity for K + would be: 


A+ = 


272. -3846.59 \fl x 10" 3 + 1.41 x 10 5 -7x 10 3 (l0- 7.191\/7xl0 3 )] 

(cm 2 /H) 


166. 


166. + 380. 


[164] 


and similarly for the OH - ion. 

Once the ionic equivalent conductivities have been determined as a function of 
the KOH concentration and temperature, the ionic diffusivities can be determined by 


(53): 


D t = 


RT\i 

F 2 


[165] 


and assuming that the Nemst-Einstein condition applies the mobilities of the ionic 
species are readily calculated from the ionic diffusivities: 


Ui 


Di_ 

RT 


[166] 


Correlations for the Solubility of Oxygen and Hydrogen Gas in KOH Solutions 

The solubilities of hydrogen and oxygen gas in KOH have been investigated by 
many investigators (80-85) where the solubility data is usually reported in the form 
of a Setschenow salt effect parameter. This salt effect parameter, k scx , can be written 
in the form 

kscx = log (cm 3 /mol) 


[167] 
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Table XVIII. Limiting equivalent ionic conductivities for different temperatures. 


A° (cm 2 /eq.-fl) 

species 

25 °C 

100 °C 

80 °C a 


OH" 

197.6 

446. 

380. 


K + 

73.5 

200. 

166. 


A 0 

271.1 

646. 

546. 



a interpolated 
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where x; is the mole fraction of gas in the electrolyte solution and x° is the mole 
fraction of the gas in pure water. For hydrogen, k scx is fairly constant at 0.129 over 
the temperature range of 25 to 100 °C . For oxygen, the salt effect parameter varies 
almost linearly over the temperature range of 60 to 100 °C such that 


k acx> o 2 = 0.1923 -0.10 x 10 3 r (cm 3 /mol) 


[168] 


where T is given in K. Bensen et al. (86) derived the following equation for the 
solubility of oxygen in water over the temperature range of 373 to 563 K: 

1.3104 x 10 4 3.4170 x 10 6 , 2.4749 x 10® 


In 


= -4.1741 + 


"02 


T 


J' 2 


+ 


r 3 


[169] 


The mole fraction of hydrogen gas in pure water was similarly fitted to experimental 
data by Battino and Wilhelm (87) over the temperature range of 273 to 353 K giving: 

5528.45 16.8893 


lnx^ 2 = —48.1611 + 


+ 


In 


[170] 


100.0 


Henry’s law constants were formulated as a function of temperature and concentration 
by rearranging Eq. [167] to give: 

= (2 C. + C') x . (t = 02i H 2 ) [171] 

1 — Xi 

where the mole fraction for species i is given by Eqs. [167] and [169] for oxygen 

1 

.1 4170x10® . 2.4749x10® 1 

[172] 


x 0 7 = 


exp 


4.1741 + llliplSl - Ml™*30i + ui 

1 

X 2 Q( 0 . 1923-0. 10xlO- 3 T)1000C e 

and by Eqs. [167] and [170] for hydrogen: 

exp [—48.1611 + 552 ^' 4 — + 16.8893 In joo] 

x n 2 = JQ129.0C,. “ 


[173] 



